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We present a complete update of the analysis of v e and v e disappearance experiments in terms 
of neutrino oscillations in the framework of 3+1 neutrino mixing, taking into account the Gallium 
anomaly, the reactor anomaly, solar neutrino data and v e C scattering data. We discuss the impli- 
cations of a recent 71 Ga( 3 He, 3 H) 71 Ge measurement which give information on the neutrino cross 
section in Gallium experiments. We discuss the solar bound on active-sterile mixing and present our 
numerical results. We discuss the connection between the results of the fit of neutrino oscillation 
data and the heavy neutrino mass effects in /3-decay experiments (considering new Mainz data) and 
neutrinoless double-/? decay experiments (considering the recent EXO results). 

PACS numbers: 14.60.Pq, 14.60.Lm, 14.60.St 



I. INTRODUCTION 

In recent years several short-baseline neutrino oscilla- 
tion experiments have found anomalies which may re- 
quire an extension of the standard three-neutrino mix- 
ing framework which describes the neutrino oscillations 
observed in solar, atmospheric and long-baseline experi- 
ments (see [lHSI)- In this paper we consider the Gallium 
anomaly [3-[6| and the reactor anomaly @~E1: which in- 
dicate that electron neutrino and antineutrinos may dis- 
appear at short distance^. Such disappearance may be 
explained by the presence of at least one massive neutrino 
at the eV scale, which drives short-baseline neutrino os- 
cillations generated by a squared-mass difference which is 
much larger than the squared-mass difference operating 
in the solar, atmospheric and long-baseline neutrino os- 
cillation experiments. We consider 3+1 neutrino mixing, 
which is the minimal extension of three-neutrino mixing 
which can explain the Gallium and reactor anomalies. 
Since from the LEP measurement of the invisible width 
of the Z boson [l3[ we know that there are only three 
light active flavor neutrinos the additional neutrino in 
the 3+1 framework is sterile. 

In this paper we discuss the implications of the recent 
71 Ga( 3 Hc, 3 H) 71 Ge measurement in Ref. 14|, which give 
information on the neutrino cross section in Gallium ex- 
periments. We take also into account the most updated 



The inclusion in the analysis of the more controversial LSND 

[Toll and MiniBooNE [TJ ( ^ -> l vl 
elsewhere 



calculation of the reactor neutrino fluxes presented in the 
recent White Paper on light sterile neutrinos [15j . We 
present also a detailed discussion of the connection be- 
tween the results of the fit of neutrino oscillation data 
and the results of /3-decay experiments (considering the 
Mainz data presented very recently in Ref. flU ]) and neu- 
trinoless double-/? decay experiments (considering the re- 
cent EXO bound in Ref. [17| and the controversial posi- 
tive result in Ref. P3). 

We consider 3+1 neutrino mixing as an extension of 
standard three-neutrino mixing. The mixing of the three 
active flavor neutrino fields v e , v^, v T and one sterile 
neutrino field v s is given by 



kVk 



(1) 



where U is the unitary 4x4 mixing matrix (U 1 = U^ 1 ) 
and each of the four 's is a massive neutrino field with 
mass Wfc. We consider the squared-mass hierarchy 



Ara^ < Awij! < At 



(2) 



i 2 , such that A7712! generates 



anomalies will be discussed 



with Amjy = m 2 . 
the very-long-baseline oscillations observed in solar neu- 
trino experiments and in the KamLAND reactor antineu- 
trino experiment, Am^ generates the long-baseline os- 
cillations observed in atmospheric neutrino experiments 
and in long-baseline accelerator and reactor neutrino and 
antineutrino experiments, and Ato^ generates short- 
baseline oscillations. 

The effective survival probability at a distance L of 
electron neutrinos and antineutrinos with energy E in 



2 



short- baseline neutrino oscillation experiments is given 
by (see Refs. [Hill) 



500 kcV 



= 1 - sin 2 2ti e 



sin 



4E 



with the transition amplitude 



sin 2 2tf ee = 4|t/ e4 | 2 



\U e 



(3) 



(4) 



The plan of the paper is as follows. In Section |TT] we 
discuss in detail the Gallium v e anomaly and the 

implications of the important recent 71 Ga( 3 He, 3 H) 71 Ge 



measurement in Ref. [14j . In Section IIIII we present 
the results of the combined analysis of Gallium data 
with reactor v e data, taking into account the reactor v e 
anomaly 0"S EBl • In Section IIVI we discuss the solar 
neutrino constraint on short-baseline v e disappearance 
p3l - l26j . In Section [V] we present the results of the global 
fit of v e and v e disappearance data, which includes also 



-N 



g.s. 



scattering data |2J7j,l28|. We con- 



front these results with the bounds on the heavy neutrino 
mass given by the data of /3-decay experiments fliil . |29| 
and neutrinoless double- j3 decay experiments [l7lll8|. Fi- 
nally, in Section IVT1 we draw our conclusions. 



II. GALLIUM ANOMALY 



The GALLEX [30H32J and SAGE [3J436J Gallium so- 
lar neutrino experiments have been tested with intense 
artificial 51 Cr and 37 Ar radioactive sources which pro- 
duce electron neutrinos through electron capture with 
the energies and branching ratios given in Tab. HI In each 
of these experiments the source was placed near the cen- 
ter of the approximately cylindrical detector and electron 
neutrinos have been detected with the solar neutrino de- 
tection reaction 



71 



Ga 



71 



Ge + e: 



(5) 



The average neutrino travelling distances are 
(i)cALLEX = 1-9 m and (L)sage = 0.6 m. The 
first line in Tab. HIl shows the ratios Rb of measured and 
expected 71 Ge event rates reported by the experimental 
collaborations. The index B indicates that the expected 
event rates have been calculated using the Bahcall cross 
sections 133 



CT B ( 51 Cr) = 58.1 x 10~ 46 cm 2 
cr B ( 37 Ar) = 70.0 x 10~ 46 cm 2 



(0) 
(7) 



TABLE I. Energy (E) and branching ratio (B.R.) of the 
neutrino lines produced in the electron-capture decay of 51 Cr 
and 37 Ar. 







51 Cr 




37 Ar 


E [keV] 


747 


752 427 


432 


811 813 


B.R. 


0.8163 


0.0849 0.0895 


0.0093 


0.902 0.098 




FIG. 1. 



Ga 



L Ge transitions induced by 51 Cr and 37 Ar 



electron neutrinos. 



without considering their uncertainties. One can see that 
the values of i?^ 1 and R^ 1 indicate a compatibility be- 
tween the measured and expected event rates, whereas 
the values of i?g 2 and R B 2 are significantly smaller than 
one, indicating a disappearance of electron neutrinos. 
The weighted average in Tab. HIl gives a 2.7<r anomaly. 

Since the values of the cross sections of 51 Cr and 37 Ar 
electron neutrinos and their uncertainties are crucial for 
the interpretation of the Gallium data as indication of 
short-baseline v e disappearance, in the following we dis- 
cuss in detail the problem of the determination of the 
cross sections and their uncertainties, taking into account 
Refs. [37| - |40| and the important recent measurement in 
Ref. [lj 

The cross sections of the interaction process ([5]) for 
neutrinos produced by 51 Cr and 37 Ar sources are given 
by 



BGT 



175 



175 



BGT„ 
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BGT 



.-.mi 



500- 



BGT„ 



(8) 



where a gs is the cross sections of the transitions from 
the ground state of 71 Ga to the ground state of 71 Ge, 
BGTg S is the corresponding Gamow- Teller strength, and 
BGT175 and BGT500 are the Gamow- Teller strengths of 
the transitions from the ground state of 71 Ga to the 
two excited states of 71 Ge at about 175 keV and 500 
keV (see Fig. EE}. The coefficients of BGTi 75 /BGT gs 
and BGT5oo/BGT gs are determined by phase space: 
W 51 Cr) = 0.669, £ 5 oo(!iCr) = 0.220, £i75( 37 Ar) = 
0-695, 6oo( 37 Ar) = 0.263 |37[. 



TABLE II. Ratios of measured and expected 71 Ge event rates 
in the four radioactive source experiments. Gl and G2 denote 
the two GALLEX experiments with 51 Cr sources [HQ, SI 
denotes the SAGE experiment with a Cr source, 
denotes the SAGE experiment with a 37 Ar source 
AVE denotes the weighted average. 



and S2 

MM- 



Rb 
Rhk 
Rff 
Rhf 



Gl 
0.95j 
0.85:1 



-+0.11 

-0.11 
■+0.12 
-0.12 

0.93±°;£ 

,+0.13 
-0.11 



G2 

o.8il° j? 

0.7lj 



0.83] 



1+0.11 
-0.11 

0.7912;!? 

1+0.11 
-0.11 



0.71] 



SI 

0.951°- 
0.841°;" 
93+ 11 

u - yo -0.12 
+0.13 
0.12 



S2 



0.83 



n 7Q+0 08 
u - ' a -0.08 
+0.09 
-0.09 

77+ - 09 

-0.10 
-0.09 



0.71] 



0.69] 



AVE 
86+ 05 

u -°°-0.05 

n 77 +o.os 

U.f '-0.08 

n s/i+ - 05 

U - 84 -0.05 
n yc+0.09 
u - ' d -0.07 
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The cross sections of the transitions from the ground 
state of 71 Ga to the ground state of 71 Ge have been cal- 
culated accurately by Bahcall [37^ : 



Ogs^Cr) = 55.3 x 10 
cr gs ( 37 Ar) = 66.2 x 10 



- 46 cm 2 



- 46 cm 2 



(9) 
(10) 



These cross sections are proportional to the characteristic 
neutrino absorption cross section [37], El[ 



cr = 2aZ Gc mlGl\V ud \ 2 g\BGT gs 

= Z Gc BGT gs (3.091 ± 0.012) x 10" 46 cm 2 



(11) 



where a is the fine-structure constant, Zq c = 32 is the 
atomic number of the final nucleus, m e is the electron 
mass, Gf is the Fermi constant, V u d is the ud element of 
the quark mixing matrix V, and <?a is the axial coupling 
constant. The numerical value of the coefficient in the 
last line of Eq. (fTTT) has been obtained with the values 
of these quantities given in the last Review of Particle 
Physics [42(. From the value 



o-o = (8.611 ±0.011) x 10 



- 46 cm 2 . 



(12) 



calculated by Bahcall [37J using the accurate measure- 
ment H3 

Ti/ 2 ( n Ge) = 11.43 ± 0.03d (13) 

of the lifetime of 71 Ge (which decays through the 
electron-capture process e~ + HGe — > l\G& + u e , which 
is the inverse of the u e detection process ©), we obtain 



BGTg S = 0.0871 ± 0.0004 . (14) 

This value agrees with that given in Ref . [401 . but it is 
different from that recommended in Ref. [l4[ . Hence, we 
checked it using the relation 



BGT„ S = 



[2J, 



Gc 



11 



2tt 3 In 2 



[2J Ga + 1] G| | V ud \ 2 mlg\ft 1 / 2 ( 71 Ge) 
6289 ± 3s 



2 5 2 A /i 1/2 ( 71 Ge) ' 
with Jq c = 1/2 and J G a = 3/2, and the value 
log /£i/ 2 ( 71 Ge) = 4.3493 ± 0.0015 , 



(15) 



(16) 



obtained with the LOGFT calculator [44[ of the National 
Nuclear Data Center using the lifetime (fT3|) . The result, 



BGT„ 



0.0872 ±0.0005. 



(17) 



is in agreement with the value (|14[) , which will be used 
in the following. 

The GT strengths BGT175 and BGT 500 have been 
measured in 1985 in the (p, n) experiment of Krofcheck 
ct al. [11 [39[ and recently, in 2011, in the ( 3 He, 3 H) ex- 
periment of Frekers et al. [l4| . The results are listed in 



Tab. IIIII together with the 1998 shell-model calculation 
of BGT 175 of Haxton 0. 

The Bahcall cross sections (j6]) and (0 have been ob- 
tained using for BGT500 the Krofcheck et al. measure- 
ment and for BGT175 half of the Krofcheck et al. upper 
limit [37 1 . 

we used the Hax- 



In previous publications [6l. l28l . l45l - l4g 
ton shell-model value of BGT175 and the (p, n) measured 
value of BGT500- Although the uncertainties of the Hax- 
ton shell- model value of BGT175 are so large that BGT175 
may be negligibly small, the central value is much larger 
than the upper limit obtained in the (p, n) experiment. 
According to Ref. p0| . this is due to a suppression of 
the (p, n) value caused by a destructive interference be- 
tween the spin (A J = 1, AL = 0) matrix element and 
an additional spin-tensor (A J = 1, AL = 2) matrix el- 
ement which operates only in (p, n) transitions. We do 
not know if the same suppression is operating also in 
( 3 He, 3 H), which could be the explanation of the small- 
ness of the value of BGT175 measured by Frekers et al., 
which is compatible with the (p, n) upper bound. More- 
over, the value of BGT500 measured by Frekers et al. has 
a 2.1a discrepancy with that measured by Krofcheck et 
al. Since we cannot solve these problems, we consider the 
following three approaches which give the cross sections 
in Tab. lIVI and the ratios of measured and expected 71 Ge 
event rates in Tab. [Til 

HK: Haxton BGT175 value and Krofcheck et al. BGT500 
value. This is our old a ppr oach adopted in previous 
publications [f| [H, l45l448| . The cross sections are 
significantly larger than the Bahcall cross sections 
in Eqs. ^ and ([7]), albeit with large uncertainties 
which make them compatible at the la level. 

FF: Frekers et al. values of both BGT175 and BGT 500 . 
This is a new approach which is motivated by the 
new ( 3 He, 3 H) measurements [l4j|. The cross sec- 
tions are only slightly larger than the Bahcall cross 
sections in Eqs. ([5]) and ©, mainly because of the 
larger BGT 500 . 

HF: Haxton BGT175 value and Frekers et al. BGT 50 o 
value. This is a new approach which is motivated 
by the possibility that the BGT175 measured by 
Frekers et al. suffers of destructive interference 
between the spin and spin-tensor matrix element 
and its value is different from the BGT175 in Gal- 
lium neutrino detection, as discussed by Haxton for 
the (p, n) experiment [401 ] - This approach gives the 
largest cross sections, which however are still com- 
patible with the Bahcall cross sections at the la 
level. 

From the weighted averages of measured and expected 
71 Ge event rates in Tab. [TTJit follows that the statistical 
significance of the Gallium anomaly in the three cases 
is, respectively, about 3. Oct, 2.9a and 3.1cr. Hence, the 
new ( 3 Hc, 3 H) cross section measurement of Frekers et 
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TABLE III. Values of the Gamow- Teller strengths of the transitions from the ground state of Ga to the two excited states 
of 71 Ge at 175 keV and 500 keV and their relative values with respect to the Gamow- Teller strength of the transitions to the 
ground state of 71 Ge, given in Eq. (|14l) . 



Reference Method 

Krofcheck et al. [38,39] 71 Ga(p, n) 71 Ge 

Haxton [40J Shell Model 

Frekers et al. [14] 71 Ga( 3 He, 3 H) 71 Ge 



R^rp BGTi 75 BGT500 

^Ol 175 T, prr BUI 500 RrT 

< 0.005 < 0.057 0.011 ±0.002 0.126 ±0.023 
0.017 ±0.015 0.19 ±0.18 

0.0034 ± 0.0026 0.039 ± 0.030 0.0176 ± 0.0014 0.202 ± 0.016 



al. [1J| confirm that there is a Gallium anomaly at a 
level of about 3c @, which indicates a short-baseline 
disappearance of v e which can be explained by neutrino 
oscillations. 

We analyzed the Gallium data in the three cases above 
in terms of neutrino oscillations in the 3±1 framework, 
in which the effective probability of v e survival is given 
by Eq. ([3]) with a = e. We used the statistical method 
discussed in Ref. [fjj], neglecting for simplicity the small 
difference between the 51 Cr and 37 Ar cross section ra- 
tios in Tab. IIVI The results of the fits are presented in 
Tab. [V] and Fig. [3] One can see that in any case neu- 
trino oscillations give an acceptable fit of the data. In 
the FF case the goodness-of-fit is smaller than in the HK 
and HF cases, because of the much smaller uncertainty 
of BGT175. The three cases give approximately the same 
best- fit value and allowed range of Am^. Instead, they 
differ in the best-fit value and allowed range of sin 2 2i9 ee : 
the FF case is in favor of smaller values of sin 2 2i9 ee than 
the HK and HF cases. 



III. FIT OF GALLIUM AND REACTOR DATA 

The reactor antineutrino anomaly_[9| stems from a new 
evaluation of the reactor D e flux 7, 8] which implies that 
the event rate measured by several reactor v e experiments 
at distances from the reactor core between about 10 and 
100 meters is smaller than that obtained without v e dis- 
appearance. This is illustrated in Fig. [3[ where we plot- 
ted the ratio R of the observed v e event rate and that 
expected in absence of v e disappearance for the Bugey- 
3J49|, Bugey-4 0, ROVN091 [U, Gosgen ILL 
[53| and Krasnoyarsk reactor antineutrino experi- 
ments. We used the reactor neutrino fluxes presented 



in the recent White Paper on light sterile neutrinos [15| , 
which updates Refs. 0-lSl- From Fig. [3] one can see that 
the reactor antineutrino anomaly has a significance of 
about 2.8cr. In the fit of reactor data, besides the above- 
mentioned rates, we consider also the 40 m/15 m spectral 
ratio measured in the Bugey-3 experiment (49| . 

The results of the fit of the reactor antineutrino data 
are presented in Tab. I VII and Fig. @] One can see that 
the preferred range of Am^ has a large overlap with 
that indicated by the Gallium anomaly, but there is a 
strong upper bound for sin 2 2"D ee of about 0.3. Therefore, 
the large-sin 2 2-d ee part of the Gallium-allowed region in 
each of the three cases considered in Fig. [5] is excluded 
by reactor data. 

The results of the combined fit of Gallium and reactor 
data are presented in Tab. lVII and Fig.0 From Fig.[5]onc 
can see that the allowed regions in the sin 2 2'i? ee -Am 2 1 
in the three cases that we have considered for the fit of 
Gallium data are quite similar, and the best-fit values of 
sin 2 2i} ee and Am^ are equal (see Tab. IVI[). This is due 
to a dominance of reactor data, which are more numerous 
and have smaller uncertainties. Hence, in the following 
we consider only the FF case, which is the one which 
is more compatible with reactor data, because it agrees 
more than the HK and HF cases with the reactor exclu- 
sion of large values of sin 2 2-)9 ee . This better agreement 
is quantified by the larger parameter goodness-of-fit (55j 
in Tab. ED 



IV. SOLAR NEUTRINO CONSTRAINT 

In this section we discuss the upper bound for sin 2 2i? ee 
which can be obtained from the data of solar neutrino 
experiments [3^, l56l - [6(| and from the data of the Kam- 



TABLE IV. Gallium cross section (in units of 10 -46 cm 2 ) and 
its ratio with the corresponding Bahcall cross section (Eqs. ([6]) 
and 0) for 51 Cr and 37 Ar neutrinos in the three cases dis- 
cussed in the text. 



TABLE V. Values of X 2 , goodness-of-fit (GoF) for 2 degrees 
of freedom and best-fit values of the 3±1 oscillation parame- 
ters obtained from the three fits of Gallium data described in 
the text. 





51 


Cr 




37 Ar 




HK 


FF 


HF 




a 


(t/ctb 


a 


(T/<7B 


Xmin 


4.8 


7.9 


4.6 


HK 


63.9 ±6.5 


1.10 ±0.11 


77.2 ± 8.1 


1.10 ±0.12 


GoF 


9.1% 


1.9% 


9.9% 


FF 


59.2 ± 1.1 


1.02 ±0.02 


71.5 ± 1.4 


1.02 ±0.02 


Amii[eV 2 ] 


2.24 


2.1 


2.24 


HF 


64.9 ±6.5 


1.12 ±0.11 


78.5 ±8.1 


1.12 ±0.12 


sin 2 2$ ee 


0.50 


0.30 


0.52 
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FIG. 2. Allowed regions in the sin 2 2i9 ee -Am| 1 plane and marginal Ax 2 's for sin 2 2$ 
fit of the results of the Gallium radioactive source experiments in the three cases discussed in the text 
corresponding to Xmm is indicated by a cross. 



and Am 2 ^ obtained from the combined 



The best-fit point 



IE 



IE 



IE 



IE 



Reactor Rates 



Bugey3-15 - 
Bugey3-40 - 
Bugey3-95 - 
Bugey4 
ROVNO 
Gosgen-38 - 



Gosgen-45 
Gosgen-65 
ILL 

Krasno-33 
Krasno-92 
Krasno-57 



Average Rate = 



20 



R = 0.93 ±0.024 



40 



60 



80 



100 



L [m] 



REA 

68.27% CL (1c) 
90.00% CL 
95.45% CL (2c) 
99.00% CL 
99.73% CL (3c) 




FIG. 3. Ratio R of the observed u e event rate and that ex- 
pected in absence of Ve. disappearance in reactor neutrino ex- 
periments. The horizontal band represents the average value 
of R with la uncertainties. 



FIG. 4. Allowed regions in the sin 2 2$ ee -Am 2 1 plane and 
marginal A\ 2, s for sin 2 2i? ee and Am 2 ! obtained from the 
combined fit of reactor antineutrino data. The best-fit point 
corresponding to Xmm i s indicated by a cross. 



LAND very-long-basclinc reactor antineutrino experi- 
ment [67} , which are sensitive to oscillations generated by 
the small squared- mass difference Am^. Since the event 
rates measured in these experiments arc well described 
by standard three-neutrino mixing, the data allow us to 
constrain the corrections due to active-sterile neutrino 
mixing, which affects the electron neutrino and antineu- 
trino survival probability and generates transitions into 
sterile neutrinos [23l - l26l . l68l [69j . As explained in the fol- 



lowing, the almost degenerate effects in the solar and 
KamLAND experiments HJ-dfj of |t/ e4 | 2 , which deter- 
mines sin 2 2i7 ee through Eq. (@|, and |t/ e 3| 2 can be re- 
solved by using the recent determination of |t/ e 3| 2 in the 
Daya Bay [7(j and RENO [7l| long-baseline reactor an- 
tineutrino experiments. 

The effective survival probability of electron neutrinos 
and antineutrinos in the solar and KamLAND cxperi- 
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10 2 10 1 2 4 6 8 10 10 2 10~ 1 2 4 6 8 10 10 2 10~ 1 2 4 6 8 10 

s/n 2 2fl 00 Ax 2 sin 2 2S„ &% 2 sin 2 2S m A% 2 



FIG. 5. Allowed regions in the sin 2 2i9 ee -Am 2 1 plane and marginal A^ 2 's for sin 2 2i} ee and Am 2 i obtained from the combined 
fit of Gallium and reactor data in the three cases discussed in Section HU The best-fit point corresponding to Xmin is indicated 
by a cross. 



merits is given bj0 [72| 



(18) 



fc=3 



k=3 



where P 2v ^ v is the two-neutrino survival probability 
generated by Am^i. Considering small values of \U e3 \ 2 
and |C/ e 4| 2 , we have 



(19) 



fc=3 



In vacuum the two-neutrino survival probability -P 2ly _>.„ 
has the standard two-neutrino form which does not de- 
pend on \U e z\ 2 and |/7 e 4| 2 - Therefore, in the KamLAND 



TABLE VI. Values of x , number of de erees of freedom 
(NDF), goodness-of-fit (GoF) and best-fit values of the 3+1 
oscillation parameters obtained from the fit of reactor an- 
tineutrino data (first column) and from the combined fit of 
reactor and Gallium data in the three cases discussed in Sec- 
tion [Il|The last three lines give the parameter goodness-of-fit 
(PG)gl of the combined fit. 





REA 


REA+HK 


REA+FF 


REA+HF 


Xmin 


21.5 


30.6 


31.8 


31.0 


NDF 


36 


40 


40 


40 


GoF 


97% 


86% 


82% 


85% 


AmlifeV 2 ] 


1.9 


1.95 


1.95 


1.95 


sin 2 2l? ee 


0.13 


0.16 


0.16 


0.16 


Axpg 




4.3 


2.4 


4.8 


NDFpg 




2 


2 


2 


GoFpg 




12% 


30% 


w 



In this discussion wc neglect, for simplicity, the matter effects in 
the Earth which affect the neutrino detection rates in the night, 
but these effects are taken into account in our calculation. 



experiment |t/e3| 2 an d |t^e4| 2 have the same effect of sup- 
pressing the electron neutrino and antineutrino survival 
probability. For solar neutrinos, the main effect of |£/ e 3| 2 
and |J7 e4 | 2 is the same as in the vacuum case, but there 
are corrections caused by the modifications of P 2 "^ Vc due 
to the decrease of \U el \ 2 + \U e2 \ 2 = 1 - {\U e3 \ 2 + |£/ e4 | 2 ) 
if | C/ e 3 1 2 7^ and/or \U e 4\ 2 ^ and to the contribution 
of the neutral current potential Vnc which is not felt by 
the sterile neutrino during propagation in matter. 

In order to describe this effect, we neglect possible CP- 
violating phases in the mixing matrix and we parameter- 
ize it as (see also [24|) 



U — R23R24R34R14R13R12 



(20) 



where R a b is the real orthogonal matrix (R^ b = R a b) 
which operates a rotation in the a-b plane by an angle 

$ab- 



[R ab ]rs = S rs + (c a b - 1) (Sra S sa + S r b S s b) 
+ S a b (Sra S s b — 5 r b Ssa) , 



(21) 



with c a b = cos i)ab and s a b = sm-& a b- In this parame- 
terization the electron line of the mixing matrix is given 
by 

U e i = c 12 c 13 c 14 , U e2 = S 12 Ci 3 C 14 , (22) 

Ue3 = S13C14 , U e 4 = S U . (23) 

Hence, this is an extension of the standard three-neutrino 
mixing parameterization of the electron line (see [H-0]) 
with the addition of v c -v± mixing parameterized by $14. 
This parameterization is also convenient because 



0i 



(24) 



Since the sterile line of the mixing matrix is more com- 
plicated, it is convenient to write its first two elements 
as 



U s \ = cos ip s cos Xs, U s2 = sin ip s cos Xs 



(25) 
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with 



tan ip s 



cos 2 Xs 



ZS12C24 + C12S24 
ZC12C24 - S12S24 



4 

E 

fc=3 



1- > \U sk \ 2 = Z 2 c 2 



24 ' ''24 



^ = C13C34S14 - S13S34 . 



(26) 

(27) 
(28) 



The adiabatic two-neutrino survival probability 
P2% is given by H 



P 2 



1 



(l + cos 2i?i2 cos 2#'j l 2 ) 



(29) 



where i?]*2 is the effective mixing angle at neutrino pro- 
duction, which is given by 



#12 + W° , 



(30) 



The mixing angle uj between the vacuum mass basis and 
the effective mass basis in matter is given by 



tan 2co 



2£V r sin2£ 



Am^ - 2£Vcos2£ 



(31) 



Here V is the matter potential given by 

^ 2 ^Cc44 + ^N C COS 4 X, 

- 2V"cc 1 1 / nc cos 2(^12 - </? s )c 2 3 c 2 4 cos 2 x 



(32) 



where Vqc and Vnc are the standard charged-current and 
neutral-current matter potentials. The angle £ is given 
by 



tan2£ 



Vcc sin27?i2C 2 3 c 2 4 - Vnc sin 2^ cos 2 Xs 
V C c cos 2i?i 2 c 2 3 c 2 4 - V^c cos 2<,c s cos 2 Xs 



(33) 

Therefore, for S14 <C 1 the contributions of | C/«=3 1 2 ~ s 2 3 
and \ U e 4\ 2 = s 2 4 to the matter effects are almost degener- 
ate. There is only a small difference of their contributions 
due to Z in Eq. pgj) . 

In solar neutrino measurements the degeneracy of the 
effects of \U e 3\ 2 and \U e 4\ 2 is also slightly broken by the 
SNO neutral current measurement, which is sensitive to 
the total probability of v e transitions into active neutri- 
nos, which by unitarity is given by 1 — P„^]^ v , with 



pSUN _ p2v 



C? 3 C 2 4 COS 2 X s + Y,\ U ' 



ek\ \U s k 



(34) 



k=3 



Here, the adiabatic two-neutrino transition probability 
P 2 % Va is given bj0 [13 

PlU„, = \ (1 + cos 2<p s cos 2tf° 2 ) . (35) 



One can check that in the limit of two-neutrino v e -u s mixing 
the unitarity relation P^ u ^ v + f >„ s = 1 is satisfied. In this 
case, C13 = C14 = S24 = cos Xs = 1 and ip s = #12 + tt/2. 



The degeneracy of |J7 e 3| 2 and |/7 e 4| 2 is broken by their 
different effects in ip s , Xs and in the last term of Eq. (|34[) . 

Luckily, the recent determination of the value of $13 in 
the Daya Bay [7(| and RENO 71 experiment allow us 
to obtain information on the value of $ e e = #14 without 
much uncertainty due to 1^13. In the 3+1 mixing scheme 
under consideration the effective long-baseline v e survival 
probability in the Daya Bay and RENO far detectors is 
given by [73j 



P, 



LBL-F 



1 



sin 2 2i?i3 sin 2 



AE 



- i sin 2 2tf 14 , 



(36) 

since the oscillations due to Ato 2 ^ 3> Am 2 ^ are averaged 
and the oscillations due to Am^ <C Ato^ are negligibly 
small. This survival probability depends on $14, but the 
Daya Bay and RENO collaborations measured the ratio 
of the probability (|3"6"|) in the far detectors and the prob- 
ability (|3T5)) with Am§ 1 L/4_E <C 1 in the near detectors, 



P- L 



LBL-N 



= 1 - 



'2t? 



14 ■ 



(37) 



Since the contribution of small values of $14 to the mea- 
sured ratio is of order i?f 4 (74|; i n practice the value of 
$13 determined by the Daya Bay and RENO collabora- 
tions with a three-neutrino mixing survival probability 
is accurate also in the 3+1 scheme under consideration. 
Nevertheless, since we have the possibility, in our anal- 
ysis we took into account the exact ratio of the far and 
near survival probabilities (|3"6")l and (f3"T|) by including the 
least-squares function Xlbl °f the far/near relative mea- 
surements of Daya Bay and RENO in the total solar and 
reactor least-squares function 



22 22 
X = Xsol + Xkl + Xlbl 



(38) 



For the calculation of the solar least-squares function 

t y , w e considered the radiochemical 37 C1 (5(| and 71 Ga 
|57| experiments, the day and night energy spectra 
of all four phases of Super-Kamiokande [58rl6l| , the day 
and night energy spectra of the SNO D2O phase (62|, the 
charged-current and neutral-current rates of the SNO salt 
63 1 and NCD 1641 phases, and the rates of low energy 7 Be 
65 1 and pep [66| solar neutrinos from the Borexino ex- 
periment. We did not use the SNO data obtained with 
the low energy threshold analysis (75| and those obtained 
with the combined analysis [76[ because both analyses as- 
sumed a three-parameter polynomial survival probability 
which is not appropriate for the sterile neutrino analysis. 
The solar neutrino fluxes are taken from the BP2004 ffl\ 
standard solar model, except for the normalization of the 
solar 8 B neutrino flux, which is considered as a free pa- 
rameter determined by the minimization of x 2 (as usual, 
because of its large theoretical uncertainties). 

The KamLAND least-squares function Xkl has been 
calculated using the energy spectrum reported in [67| 
with a total exposure of 3.49 x 10 32 target-proton-year. 

In our analysis, we used the v e survival probability (118[) 
and the v e — > v a transition probability (|34l) taking into 
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FIG. 6. Marginalized Ax 2 = X 2 ~ Xmin as a function of 
sin 2 2i9 ee obtained from the fit of solar and KamLAND data 
with and without Daya Bay and RENO data. 



FIG. 7. Allowed 95% CL regions in the sin 2 2tf ee -Aml 1 plane 
obtained from the separate fits of Gallium, reactor, solar and 
v e C scattering data and from the combined fit of all data. The 
best-fit points corresponding to Xmin are indicated by crosses. 



account as parameters the squared-mass difference Am^ 
and the five relevant mixing angles $12, #13, #14, $24, 
$34 (solar and KamLAND oscillations are independent 
from $23 j because and v T are indistinguishable; see 
[l[). The six-dimensional parameter space is explored 
with a Markov Chain Monte Carlo sampling in order to 
minimize the total \ 2 m Eq. (|38|) . 

Since the best fit is obtained for #14 = 0, the marginal- 
ized Ax 2 = X 2 ~ Xmin shown in Fig. |6]give stringent con- 
straints on the value of sin 2 2i} ee = sin 2 2i?i4. In Fig. [5] 
we have plotted the Ax 2 obtained with and without in- 
cluding the Daya Bay and RENO data. One can see that 
these data are useful in order to tighten the upper bound 
on sin 2 2$ ee . 



V. GLOBAL FIT 

In this section we present the results of the global fit of 
electron neutrino and antineutrino disappearance data, 
which includes the Gallium and reactor data discussed 
respectively in Sections [TT] and IIII[ the solar neutrino 
constraint discussed in Section IIV1 and the KARMEN 
[Zl |z| and LSND gtj v e + 12 C -4 12 N g . s . + e" scatter- 
ing data [I?) , with the method discussed in Ref. [28| . 

Figure [7] shows a comparison of the allowed 95% CL 
regions in the sin 2 2$ ee -Am 2 :1 plane obtained from the 
separate fits of Gallium, reactor, solar and v e C scattering 
data and from the combined fit of all data. One can see 
that the separate allowed regions overlap in a band delim- 



ited by Am 2 41 > 1 cV 2 and 0.07 < sin 2 2i} ee < 0.09, which 
is included in the globally allowed 95% CL region. Fig- 
ure [5] shows the globally allowed regions in the sin 2 2i9 ee - 
Am 2 ! plane and the marginal Ax 2 's for the two oscilla- 
tion parameters. The best-fit point is at a relatively large 
value of Ato 2 1; 

(Am 2 1 )bf = 7.6 eV 2 , (sin 2 20 ee )w - 0.12 , (39) 

with xj^/NDF = 45.5/51, corresponding to a 69% 
goodncss-of-fit. However, there is a region allowed at 
la around Ara\ Y ~ 2cV 2 and sin 2 2d ee ~ 0.1. The slight 
preference of the global fit for Am^ ~ 7.6 eV 2 with re- 
spect to Ato^ ~ 2 eV 2 (see the marginal A% 2 for Am^ 
in Fig. [8]), which is preferred by Gallium and reactor data 
(see Tabs. El and EH and Figs. [2j Hand EJ), is due to the 
u e C scattering data, whichprefer larger values of Amfj 
(see the discussion in Ref. |2a|). 

Comparing the minimum of the % 2 of the global fit 
with the sum of the minima of the \ 2 of the separate fits 
of Gallium, reactor, solar and v e C scattering data, we 
obtained Axp G = 11.5, with 5 degrees of freedom, which 
gives a parameter goodness-of-fit of 4%. Therefore, the 
compatibility of the four data sets is acceptable. 

The results of the global fit, as well as the results of the 
fits of Gallium and reactor data, lead to lower limits for 
Am 2 !, but there is no upper limit for Ara^ in Figs. 
IU [5] and [H Hence, one can ask if there are other mea- 
surements which constrain large values of Am^. The 
answer is positive and comes from the measurements of 
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FIG. 8. Allowed regions in the sin 2 2i9 e e-Am| 1 plane and 
marginal Ax 2 's for sin 2 2i9 ee and Am 2 ! obtained from the 
global fit of i/ e and i> e data. The best-fit point corresponding 
to Xmin is indicated by a cross. 



FIG. 9. Relative deviation of the Kurie plot in /3-decay for 
some points in the allowed regions of Fig. [8] 



the effects of heavy neutrino masses on electron spec- 
trum in /3-decay far from the end-point, from the results 
of neutrinoless double-/3 decay experiments for Majorana 
neutrinos, and from cosmological measurements. In this 
discussion we consider the mass hierarchy 



7714 » 7771, 7772, 7713 , 



which implies 



7774 



(40) 



(41) 



Let us consider first /3-decay experiments. The ratio 
of the Kurie function K(T) in /3-decay for the case of 
a heavy neutrino ^4 and that corresponding to massless 
neutrinos is given by [H] 



K{T) N 
Q-T, 



1- \UeA? 



1- 



(Q - ry 



9(Q -T- 7774) , 



(42) 



where T is the kinetic energy of the electron, Q = 
18.574 keV is the Q-value of the decay 9 is the Hcaviside 
step function, and we have neglected the contribution of 
the three light neutrinos V\ , 1/2, v%. Figure [3] shows the 
relative deviation of the Kurie plot with respect to the 
massless case for some points in the allowed regions of 
Fig. [5J One can see that in order to see the effect of 
7774, beta-decay experiments must have a sensitivity to 



the relative deviation of the Kurie plot of the order of a 
percent or better for T > Q — 7774. 

In 2001 the Genoa 187 Re beta-decay experiment [29[ 
searched for deviation of the electron spectrum due to a 
heavy neutrino with a mass from 50 to 1000 eV. From 
Fig. 3 of Ref. one can see that the 95% CL up- 
per bound for 7774 is about 300 eV if sin 2 2i? ee ~ 0.1, 
which implies a very large upper limit on Amf^ of about 
10 5 eV 2 . 

Very recently, the Mainz collaboration released new 
data obtained with the phase II of the Mainz Neu- 
trino Mass Experiment [l6| which constrain the value of 
sin 2 77 ee for m 2 between about 10 and 3 x 10 4 eV 2 . Fig- 
urelTOl shows the constraints in the sin 2 27? ee ~Am| 1 plane 
that we obtained with a x 2 analysis of the Mainz data in 
[l6j |. From the comparison with the allowed regions ob- 
tained from the global fit of v e disappearance data shown 
in Fig. [10] one can see that the Mainz data constrain 
Am 2 ! to be smaller than about 10 4 eV 2 at about 90% 
CL. This is confirmed by the results of the combined fit 
shown in Fig. 1111 

The KATRIN experiment (see [8lt|). which will start in 
2015 [13], may be able to improve dramatically the upper 
limits on 7774 and maybe see its effects on the electron 
spectrum [83j. 

The heavy neutrino mass 7774 has also an effect in neu- 
trinoless double- /3 decay (see |84T [87 | n . if massive neu- 
trinos are Majorana particles (see [l|-|3(). Considering 
Eq. (|4ip . the contribution of the heavy neutrino mass 
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FIG. 10. Comparison of the Mainz /3-decay bound (curves 
in the top-right part of the figure) with the allowed regions in 
the sin 2 2i9 ee -Am4 1 plane obtained from the global fit of v e 
disappearance data (same as in Fig. [S]). 



FIG. 11. Allowed regions in the sin 2$ ee -Am± 1 plane ob- 
tained from the combined fit of u e disappearance and Mainz 
/3-decay data. 



TO4 to the effective Majorana mass 



is given by 



,( 4 ) 



\u, 



c4 



(43) 



(44) 



Figure [T2"l shows the marginal Ax 2 = X 2 ~ Xmm as a func- 
tion of rn^gl obtained from the global fit. One can see that 



,(4) 



rrigi is bounded from below and it is likely to be larger 
than about 10 -2 eV, a value which may be reached in 
the next generation of neutrinoless double-/? decay exper- 
iment (see (UHD]). Of course, if the three light neutrinos 
fx, i>2, v?, are quasi-degenerate at a mass scale larger than 
about 10 -2 eV the contribution of can cancel with 
that of the three light neutrinos. Such an unfortunate 
cancellation can also happen if the masses of the three 
light neutrinos follow an inverted hierarchy [89l . |90T | , since 
in that case their contribution to the effective Majorana 
mass is [85l 



1.4 x IO" 2 < m ( ii ght) < 5.0 x 10~ 2 eV 



(IH-95%CL). 

(45) 

On the other hand, no cancellation is possible in the case 
of a normal hierarchy, for which 



oft t] ~ 4 - 5 x 10 " 3 cV ( NH _ 95% CL ) ' ( 46 ) 



and the contribution of is dominant if it is larger 
than about 10 _2 cV. Figure [TBI shows the allowed values 
of mpp as functions of the lightest neutrino mass in the 
two 3+1 schemes with a normal (left) and inverted (right) 
mass spectra of the three light neutrinos. We have drawn 
also the curves which delimit the three-neutrino mixing 
allowed regions [85j. One can see that practically the 
situation is reversed with respect to the three-neutrino 
mixing case (see also the discussion in in which mpp 
is predicted to be large in the inverted spectrum and can 
vanish in the normal spectrum. In the 3+1 case mpp can 
have any value in the inverted spectrum, whereas in the 
normal spectrum it is likely to be large if the three light 
neutrino masses are hierarchical, i.e. if mi <C rri2 <C m^. 
Let us consider the "no-cancellation" case, in which 



nip 13 > TO /3/3 (no-cancellation) 



(47) 



In this case, as shown in Fig. 1121 large values of tti^I are 
excluded by the currently most stringent upper bound for 
mpp obtained in the EXO experiment [171 ] (the vertical 
green band in Fig.UJis the 90% CL EXO bound taking 
into account nuclear matrix element uncertainties). This 
limit implies the upper bound on Am^ ~ (m^/|[/ e 4| 2 ) 2 
as a function of sin 2 2i? ee shown in Fig. [TU One can 
see that parts of the high-Am 2 ^ regions allowed by the 
global fit are disfavored by the EXO bound. However, 
the large nuclear matrix element uncertainties do not 
allow to establish a precise bound. From Fig. [JJ] one 
can also see that the putative Klapdor-Klcingrothaus et 
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FIG. 12. Marginal A\ 2 = X 2 — Xmin as a function of m,pp 
obtained from the global fit. The vertical green band repre- 
sents the currently most stringent 90% CL upper bound for 
rrSpp in the no-cancellation case (|47[) obtained from the 90% 
CL EXO bound on map taking into account nuclear matrix 
element uncertainties [171 ] . The vertical orange band corre- 
sponds to the la Klapdor-Kleingrothaus et al. range of mpp 

r 



al. la range of nipp [l8[ implies a rather large value of 
Am];, around 100 — 200 eV 2 . In this case the oscillation 
length is very short, of the order of 1 cm for neutrinos 
with energy of the order of 1 MeV, as reactor neutrinos 
and neutrinos emitted by radioactive sources. Hence it 
will be practically impossible to observe a variation of 
the event rate characteristic of oscillations in future very 
short-baseline reactor neutrino experiments l9ll [§3 and 
radioactive source experiments [9J; 



496J (see also [15|, l9J ) . 
Finally, considering cosmological measurements one 



must say that they are a powerful probe of the number 
of neutrinos and of neutrino masses at the eV scale (see 
[FoI [§9[ ) , but the analysis requires many assumptions 
on the cosmological model and its details. A comparison 
of the results of the fit of short-baseline oscillation data 
is beyond the scope of this paper. We can only say that 
the analysis of cosmol ogical dat a in the framework of the 
standard ACDM [H [lOOl - fToBj allow the existence of a 
sterile neutrino thermalized in the early Universe, but 
restricts its mass to be less than about 1 eV (a possible 
suppression of the sterile neutrino thermalization with a 
large lept on asymmetry has been discussed recently in 
[l06lll07l |). If this constraint is correct, the upper bound 
on Am]] is about leV 2 , which is much more restrictive 
than those of /3 decay in Fig. [TU] and neutrinolcss double- 
/3 decay in Fig. [TU 

VI. CONCLUSIONS 



In this paper we presented a complete update of the 
analysis of v e and v e disappearance experiments in terms 
of neutrino oscillations in the framework of 3+1 neutrino 
mixing. We have shown that the Gallium anomaly, the 
reactor anomaly, solar neutrino data and v e C scattering 
data are compatible with short-baseline oscillations with 
an amplitude sin 2 2$ ee between about 0.03 and 0.2 and a 
squared-mass difference Am]] larger than about 0.5 eV 2 
at 95% CL. Assuming the mass hierarchy in Eq. ([4"0"|). we 
have shown that the heavy neutrino mass is observ- 
able in /3-dccay experiments and neutrinolcss double-/? 
decay experiments. The very recent Mainz /3-decay data 
[lj| constrain Am^ to be smaller than about 10 4 eV 2 at 
95% CL. For Majorana neutrinos, the recent EXO limit 
on the effective Majorana mass in neutrinoless double-/3 
decay [l7j give a more stringent constraint which can vary 
between about 10 2 and 10 3 eV 2 depending on the nuclear 
matrix clement uncertainties if there are no cancellations 
between the contribution of v\ and that of the three light 
neutrinos. We think that our results are interesting for 
the many projects which will search in the next years ef- 
fects of light sterile neutrinos with electro n neutrin o and 

iJJjMIIil) and 



antincutrino radioactive sources (see 
reactor electron antineutrinos (see [l5 
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